Designing short robust NOT gates for quantum computation 



CO 

i— H 

o 

<N 



> 

a 

<N 
> 

t> 

cn 
o 

CO 



Jonathan A. Jones 1 B 

1 Centre for Quantum Computation, Clarendon Laboratory, 
University of Oxford, Parks Road, 0X1 3PU, United Kingdom 
(Dated: May 8, 2013) 

Composite pulses, originally developed in Nuclear Magnetic Resonance (NMR), have found 
widespread use in experimental quantum information processing (QIP) to reduce the effects of 
systematic errors. Most pulses used so far have simply been adapted from existing NMR designs, 
and while techniques have been developed for designing composite pulses with arbitrary precision 
the results have been quite complicated and have found little application. Here I describe techniques 
for designing short but effective composite pulses to implement robust not gates, bringing together 
existing insights from NMR and QIP, and present some novel composite pulses. 



PACS numbers: 03.67.-a, 82.56.-b 

Quantum Information Processing (QIP) is the encod- 
ing of information in two level quantum systems called 
qubits, and the manipulation of this information through 
a series of unitary transformations which can be inter- 
preted as logic gates [1]. Building real quantum com- 
puters will require the ability to perform accurate uni- 
tary transformations on quantum systems in the presence 
of realistic errors. Such errors can be divided into two 
broad categories: random errors, arising from decoher- 
ence, which can be tackled by methods such as quantum 
error correction (HQ, and systematic errors, arising from 
imperfections in control fields. If systematic errors vary 
slowly in time, or if they vary over a spatial ensemble of 
qubits, it is necessary to design control sequences which 
are intrinsically tolerant of a range of error values. 

One particularly successful approach, adopted from 
Nuclear Magnetic Resonance (NMR) experiments, is the 
use of composite pulses (j-tZj, in which a single rotation 
about some axis in the x?/-plane is replaced by a sequence 
of such rotations, such that the combined propagator im- 
plements the desired rotation in the absence of errors, 
while in the presence of small errors the errors in indi- 
vidual rotations do not accumulate but instead mostly 
cancel one another. Many traditional composite pulses 
widely used in NMR are not suitable for QIP, as they are 
"point to point" pulses which make assumptions about 
the initial and final states. For example, composite inver- 
sion pulses (8l4Toj are designed to interconvert the com- 
putational basis states, but do not act correctly on su- 
perposition states. However so-called Class A composite 
pulses, or general rotors, which are error tolerant for any 
initial state, are suitable 0-0 • 

Here I will concentrate on attempts to construct not 
gates, that is tt rotations about the x-axis of the Bloch 
sphere, using only tt rotations around axes in the xy- 
plane. Composite pulses of this kind have many desirable 
properties, which will allow design techniques to be ex- 
plored while sidestepping many complexities. Although 



not gates find only limited application in theoretical QIP 
they play an important role in many experimental tech- 
niques, such as dynamical decoupling [lll-[l3j. As with 
other composite pulses developed in NMR these can be 
applied in a wide range of other experiments [l4]-[l8| . 
I will describe each rotation (pulse) by its propagator 



Oct, = cos(0/2)l - isin(<9/2)cr 



(1) 



(2) 



where is the rotation angle and 0, the pulse phase, fixes 
the rotation axis in the xy-plane (note that <f> is only 
defined up to multiples of 2tt). As these are propagators 
a sequence of pulses must be written with time running 
from right to left, the reverse of the usual order for pulse 
sequences. I will also use the notation z a for a rotation 
by an angle a around the z-axis. 

It is sometimes convenient to characterize such pulses 
by their propagator fidelity ^ 



F= \ti(U ] V)\/ti(U ] U) 



(3) 



where V is the propagator of the composite pulse in the 
presence of errors and U is the desired unitary transfor- 
mation; equivalently pulses can be categorized by their 
infidelity, defined by 1 = 1 — T . (Note that this fidelity 
is essentially the Hilbert-Schmidt inner product between 
U and V.) One traditional approach is to expand the 
fidelity as a Taylor series in the size of the underlying 
error term, and then to seek to set as many low order co- 
efficients to zero as possible [19] . Alternatively one may 
isolate the error term in the propagator itself and then 
expand this as a Taylor series [20l422j. 

Starting from two identities (here and elsewhere I ne- 
glect physically meaningless global phases) 



K/3 K a = Z2((3-a) z (3 = 7T a Z-fi 



7T, 



'a+0/2 



(4) 
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it can be immediately deduced that any sequence of 
n individual tt rotations corresponds either to some z- 
rotation (for even n) or to some tt^ rotation (for odd n). 



In the latter case 

n 
J = l 

and so to implement a not gate the individual pulse 
phases should be chosen such that $ = 0. 

I. FIRST ORDER ERRORS 

As a simple example consider pulse strength errors, 
which occur when the strength of the driving field used 
to induce transitions between qubit states deviates from 
its nominal value by some fraction e. In this case the 
rotation angle of each pulse is also increased by some 
fraction e so that 

V = [(1 + €)tt]0 =7r (e7r) (6) 

giving 

tfV = (7re) = 1 - e(i7r/2)^ + 0(e 2 ) (7) 

so that the single pulse contains a first order error term. 
The corresponding fidelity is given by 

T = | cos(e^/2)| = 1 - eV/8 + 0(e 4 ) (8) 

and so the pulse has a second order infidelity term. In 
general a pulse with an error term of order n will have 
infidelity of order 2n. 

In the presence of pulse strength errors the total prop- 
agator can be written as 

with S = e7r. Any composite pulse of this form will equal 
the desired propagator up to order zero as long as the 
phases are chosen such that = 0, but to consider the 
first order error term it is more convenient to isolate the 
error term from the ideal evolution. The identity 

6(3 7T a = 7T a 2a -(3 (10) 

allows the tt pulses to propagated to one end of the se- 
quence 

V = 7T 0ri 7T ( f) ri _ 1 . . . 7T0! 6p n . . . (fy/ (11) 

where the modified phases known in NMR as the 
phases in the interaction frame [4] or toggling frame [23|, 
l24j . are given by 

<P> = (-1) U+1) <P, + ^(-l) (fe+1) 20 fe . (12) 

k<j 

Alternatively these phases can be described recursively 
#+i - # = (-l) J (^+i - 0i) 0i = 01 (13) 
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FIG. 1: Vector diagrams showing the cancellation of first or- 
der pulse strength error terms in composite pulse sequences 
with (a) n = 3 and (b) n — 5. The equilateral polygons can 
be traced in either direction, depending on the precise rela- 
tionship chosen between phases in the toggling frame, but the 
overall orientation is fixed by the constraint on 

to give a form which is more useful in some cases. 
Expanding the error terms to first order 

5 a ^t-i(5/2)o- a + 0(5 2 ) (14) 

allows the combined error term to be approximated as 

3 

and so a composite pulse which corrects first order pulse 
strength errors can be found by choosing the such that 

5>*} = °- ( 16 ) 

3 

Since the operators in this sum all lie in the xy-plane 
and all have the same size, they can be mapped onto 
two-dimensional unit vectors and the sum solved geo- 
metrically: the underlying vectors must form a closed 
equilateral polygon of order n, as shown in Fig. HJ 

The geometric approach is particularly effective when 
n = 3 as the equilateral triangle is uniquely defined up 
to rotations and reflections. In particular such triangles 
require that 

03 = <j/ 2 ± 2tt/3 and <j>' 2 = </>[ ± 2tt/3 (17) 

where the two ± signs must be chosen with the same 
sense. Together with the requirement that <I> = these 
lead to three simultaneous equations which can be writ- 
ten as 

/0 -1 1\ /M /±2tt/3\ 

1; ?J (t) - b /3 J m 

with solutions 

(0) = <h, fo) = (T27T/3, t4tt/3, t2tt/3). (19) 

A second pair of solutions can be found by using an alter- 
native third constraint, namely $ = 7r; this relies on the 
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FIG. 2: (Color online) Fidelity achieved by (a) a simple pulse and composite pulses with n — 3 optimized to suppress first-order 
pulse strength errors (b) and off-resonance errors (c). Fidelity is plotted as a function of the fractional pulse strength error e 
and the off-resonance fraction /. Contours are drawn at 90%, 99% and 99.9% fidelity, that is logarithmically spaced infidelities 
with the inmost contour at an infidelity of 10 -3 . 



observation that tt 7T =7To up to irrelevant global phases. 
All four of these pulse sequences have the same fidelity 

T = 1 - 37r 4 e 4 /128 + 0(e 6 ) (20) 

and are entirely equivalent to each other. 

The second important type of systematic error in NMR 
pulses is off-resonance errors. These occur when the driv- 
ing field is not exactly in resonance with the underlying 
transition between the spin states. In this case 

V{0, <t>) = exp[-i#(a + fa z )/2] (21) 

where and now describe the behaviour of the pulse on- 
resonance and /, the off-resonance fraction, is the ratio 
between the frequency offset and the rotation frequency 
induced by the driving field. In this case the situation 
is more complex than was seen for pulse strength errors, 
as the propagator is no longer easily divided into desired 
and undesired terms. However it is possible to expand 
the propagator to first order, giving for a tt pulse 

V(tt,0) «7r A 0+7r/2 (22) 

with A = 2f. As before, the tt pulses can be propagated 
to one end of a composite pulse sequence 

V = 7T0 n 7T ( f )n _ 1 . . . 7T 01 X<f>» X^_ ± . . . X^ (23) 

where the modified phases are now given by 

$ = <P'j + tt/2. (24) 

As before a composite pulse which corrects first order off- 
resonance errors can be found by choosing the 0" such 
that 

3 



with the underlying vectors forming an equilateral poly- 
gon. Once again the case n = 3 is easily solved, to give 
01 = 03 = tt/3 and 02 = 27r/3, and three other closely 
related solutions, which have fidelity 

.F=l-p±^)/ 4 + 0(/ 6 ). (26) 

These composite pulses have been known for many 
years: see for example 0, an d n °te that the se- 
quence for correcting pulse strength errors was subse- 
quently generalized as the SCROFULOUS family of com- 
posite pulses [5j. The properties of the first order error 
correcting sequences derived so far are summarized in 
Fig.[2j For applications in conventional NMR it is usually 
desirable to have moderately effective error suppression 
at large errors, but for applications in QIP it is necessary 
to achieve very precise rotations for small or moderate 
errors. For this reason, one should concentrate on the 
area enclosed by the inmost contour line, corresponding 
to 99.9% fidelity. Note that in conventional NMR treat- 
ments the first-order error discussed here is referred to 
as the zero-order error, and this numbering convention 
continues at higher orders. 

Examining the inner contour lines it is immediately 
clear that it is possible to greatly broaden the high pre- 
cision region in one dimension, but that this is bought at 
some cost in increased sensitivity in the other dimension. 
Later I will consider the possibility of either creating even 
greater precision in one dimension or of broadening this 
region in both dimensions. Both of these aims will require 
sequences with at least five pulses. 

II. GEOMETRY 

It is instructive to consider an alternative "brute-force" 
approach for tackling this problem. Consider a general 
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sequence of three pulses with $ = 

(0) = (a, 0, P-ol). (27) 

In the presence of pulse strength errors the fidelity of 
such a composite pulse has the form 

2 2 

T = [3 + 2coso; + 2cos(a- / 3) + 2cos/3]+0(e 4 ) 

(28) 

and a good composite pulse can be found by minimis- 
ing the trigonometric term (note that the form of the 
fidelity expression guarantees that this term will be non- 
negative, and so seeking to minimize it is equivalent to 
seeking to set it to zero). Plotting it as a function of 
a and /3 indicates that the minima lie along the line 
f3 = — ce, and imposing this constraint allows the solu- 
tion to be easily located at a = ±2tt/3. 

Although this might seem quite different from the pre- 
vious approach, they are in fact very similar. The second 
order infidelity term in the Taylor series is simply given 
by half the magnitude of the first order error term in the 
propagator. Clearly setting this term to zero and setting 
its magnitude to zero achieves the same effect. How- 
ever looking directly at the individual error terms allows 
geometric insights to be used, greatly reducing the com- 
plexity of the calculations. This advantage becomes less 
clear once higher order terms are considered, and a judi- 
cious combination of geometry and algebra may be the 
best way to proceed. 

The role of geometry in error tolerance was recently 
discussed by Ichikawa et al. [26[. They were interested 
in distinguishing between Class A composite pulses and 
point-to-point pulses by interpreting Class A pulses as 
geometric quantum gates, that is gates where the dy- 
namic phase vanishes for every eigenstate of the ideal 
gate operator [26M30|. They showed [26| that for com- 
posite pulses where the first order error defined above 
vanishes the dynamic phase also vanishes, and so every 
Class A composite pulse must also be a geometric quan- 
tum gate. 

This observation only applies, however, to the suppres- 
sion of the first order error term, and it is necessary to 
consider separately how higher order errors can be sup- 
pressed. Note also that even for Class A pulses the error 
suppression will not be the same for every initial state: 
the composite pulse described above suppresses the first 
order pulse strength error for every initial state, but also 
suppresses the second order error for basis states. In 
other words a composite pulse can be better as an inver- 
sion pulse than as a general rotor. The fidelity measure, 
Eq. 03 in effect reports the fidelity for the worst case ini- 
tial states. 

The geometric approach was also adopted by Merrill 
and Brown 0], who noted the geometrical vector inter- 
pretation of the first order error terms. They then pro- 
ceed, however, to interpret higher order terms geometri- 
cally by considering the underlying Lie algebra. 



III. SYMMETRY 

The composite pulse sequences derived above are time 
symmetric, that is the sequence of pulse phases is the 
same when reversed. The fidelity plots in Fig. [2] are also 
symmetric around the line / = 0. These two facts are 
related: it can be shown that for any time symmetric 
sequence of tt pulses the fidelity is an even function of 
the off-resonance fraction /, containing only even powers 
in its Taylor expansion. 

In general, however, composite pulses need not be sym- 
metric, and in this case the fidelity function will contain 
odd powers in its Taylor expansion, and the fidelity plot 
will not be symmetric. This asymmetric plot can be re- 
versed around / = either by reversing the order of all 
the pulses, or by negating the phases of all the pulses. For 
a time antisymmetric composite pulse, with <\> n = — </>i 
and so on, these two operations are equivalent. For a 
time symmetric pulse they are different, but neither af- 
fects the symmetric fidelity function (3lj . 

The response to pulse strength errors is different. In 
the absence of off-resonance errors the fidelity of a se- 
quence of 7r pulses is always an even function of the frac- 
tional pulse strength error e. This can be seen by noting 
that a negative value of e is equivalent to a positive value 
of e combined with a shift in the pulse phase by tt. In the 
presence of off-resonance errors, however, this argument 
breaks down, and no particular symmetry is found. 

The symmetry of composite pulses also has effects in 
the toggling frame: a composite pulse with antisymmet- 
ric phases <pj always has symmetric phases <p f - , a fact that 
will be of some importance later. The converse does not 
quite apply: a composite pulse with symmetric phases 
need not have antisymmetric phases in general, but if 
the pulse has $ = 0, so that it is correct to zero order, 
then it will do so. As I only consider pulses which are cor- 
rect to zero order this relationship will be assumed from 
now on. Note that antisymmetric sequences of tt pulses 
always have & = 0, and so provide plausible candidates 
for composite pulse not gates [19[. 

For an antisymmetric pulse the central phase is fixed 
at 0( n +i)/2 = 0, while for a symmetric pulse it is appar- 
ently free; however the requirement that = fixes this 
value, and so both symmetric and antisymmetric com- 
posite pulses have (n — l)/2 free phases which can be 
adjusted. The more general non-symmetric composite 
pulses have n — 1 free phases, twice as many. However 
the simpler symmetric and antisymmetric pulses have 
many advantages and will be seen frequently in the sec- 
tions that follow. Antisymmetric composite pulses have 
particular advantages in some conventional NMR exper- 
iments [H, [33[ , where it is possible to render effectively 
invisible any signal arising from certain types of error, so 
that these errors lead to a drop in overall signal strength 
but do not generate any erroneous signal terms. This 
phenomenon is not normally useful, however, in QIP ex- 
periments, where symmetric composite pulses are often 
preferred. 
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IV. SECOND ORDER CORRECTIONS 

The composite pulses described above provide the best 
correction of pulse strength errors possible with a se- 
quence of three pulses, but it is possible to also remove 
second order errors with a sequence of five pulses. The 
most obvious approach is to expand the error propagator 
up to second order, but this is not ideal as the second 
order terms in this non-linear propagator include contri- 
butions from the first order error [24[ , and it is convenient 
to separate these from the genuinely second order terms. 

The traditional approach in the NMR community is 
to use average Hamiltonian theory QIJ [24, 34, 35], but 
for composite pulses it can be simpler to proceed directly 
from the underlying Baker-Campbell-Hausdorff (BCH) 
relation [35|, [36| 

e B e A = exp{£ + A + \ [B, A] + . . . } (29) 

where terms above second order have been dropped. 
This can be extended in the obvious way, and for pulse 
strength errors leads to 

S K S K _ x . . . 5 01 = exp{Ax + A 2 + . . . } (30) 

with 

ai=-4e^ ( 3i ) 

3 

and 

3 k<j 

= ij££sin(# (33) 

3 k<j 

We have already seen how to set Ai = 0, removing the 
first order error, while to remove the second order error 
requires that 

3 k<j 

A traditional NMR approach to achieve this is to 
choose an antisymmetric pulse sequence, with a sym- 
metric sequence of phases in the toggling frame, as this 
symmetry forces the terms in A 2 to cancel. Thus any 
antisymmetric composite pulse which suppresses the first 
order pulse strength error will automatically suppress the 
second order error. It remains to be shown that an 
antisymmetric composite pulse with no first order er- 
ror can be found in the case n = 5. In fact a unique 
pair of solutions exist, with 0i = 30, 02 = _0 an d 
= ± arccos(— 1/4), the well known Fi pulse pH 137^ 
(note that the order of labelling phases 1 and 2 is re- 
versed in some descriptions). The vectors describing the 
corresponding toggling frame phases, <\>' x = 0g = 3-0, 



02 = ^4 = and 03 = 4-0, form an equilateral pen- 
tagon, but it is certainly not a regular pentagon. In fact 
it is a stellated pentagon, with one side passing through 
a vertex, as shown in Fig. [H 

The Fi composite pulse has a fidelity 

T = 1 - 57T 6 e 6 /1024 + 0(e 8 ). (35) 

This fidelity can be equalled by a number of other com- 
posite pulses, but no five pulse sequence is significantly 
more effective for suppressing pulse strength errors. Its 
response in the presence of simultaneous pulse strength 
and off-resonance errors is shown in Fig. [3]^a). This plot 
shows the solution with a positive value of ip; the sec- 
ond solution is simply the time-reversed (negative phase) 
equivalent, and has an equivalent (but mirror image) fi- 
delity. Note that in this figure the inmost contour is at 
99.99% fidelity, ten times better than the inmost contour 
in Fig. [21 



V. REORDERING PULSES 

The phase angles in the Fi pulse are closely related 
to the better known BB1 pulse [H, [Hj], and here I con- 
sider the possibility of reordering sequences of 7T0 pulses. 
The fidelity definition, Eq. [3j is invariant under cyclic re- 
orderings of its components, and Eq. [10] allows the ideal 
7To pulse to be moved through the fidelity definition as 
long as the corresponding phases are negated [19]. 

Taken together these results lead to the conclusion that 
in the presence of pulse strength errors the fidelity of the 
Fi composite pulse 

(0> Fl ) = (30, 0, 0,-^,-3^) (36) 

is identical to the fidelity of the reordered pulse sequence 

(0bbi) = (0, 0, 30, 30>, 00- (37) 

This argument does not, however, apply to off-resonance 
errors, and the response of the two pulse sequences above 
is very different. Neglecting the initial 7To pulse the BB1 
pulse sequence comprises a nested set of 2tt rotations, 
and to first order such rotations have no off-resonance 
error. This can be seen either by expanding the single 
pulse propagator directly [19], or by using the first order 
expansion of a 7r pulse, Eq.[22j and noting that the second 
7r pulse generates a spin echo (39[, which removes the 
first order error terms. As a consequence the first order 
response of a BB1 pulse to off-resonance errors arises 
solely from the initial 7To pulse, and so the response is 
very similar to a simple pulse. The similarity can be 
increased still further by splitting the tto pulse into two 
and moving one half to the end of the sequence [19[. 

This fully time symmetric BB1 sequence, with the 
propagator 

(7r/2) 7r^(27r) 3 ^7r^(7r/2)o (38) 
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FIG. 3: (Color online) Fidelity achieved by some composite pulses designed to suppress pulse strength errors: (a) Fi, (b) 
symmetric BB1, (c) a new symmetric pulse with n — 5. The inmost contour is drawn at an infidelity of 10 -4 . 



has a frequency symmetric response to off-resonance er- 
rors, so the fidelity is purely an even function of /. The 
response of BB1 to simultaneous errors is compared with 
Fi in Fig.[3fb). Although BB1 does not actively suppress 
off-resonance errors, the increased sensitivity to such er- 
rors exhibited by Fi is overcome by symmetrisation. 



VI. SYMMETRIC SEQUENCES 

While this BB1 composite pulse is time symmetric, it 
breaks the informal rule of considering only sequences of 
7r pulses. It is interesting to consider whether an intrin- 
sically time symmetric sequence of five tt pulses exists 
which suppresses pulse strength errors to second order. 
In this case it is necessary to consider both the first order 
and second order errors, but these can be taken one at 
a time. In the discussions above Eq. [16] has been inter- 
preted geometrically, but it can instead be considered as 
the pair of equations 



(39) 



For a time symmetric sequence the toggling frame phases 
will be antisymmetric, and so the sum of sin terms above 
will automatically be equal to zero. Thus for the general 
symmetric sequence 

(</>) = {a,P,2P-2a,P,a). (40) 

the first order error will be suppressed if 

1 + 2 cos(a) + 2 cos(2a - 0) = 0. (41) 
This has the solution 

P = 2a±arccos[-(l + 2cosa)/2] (42) 

which is only defined when 7r/3<a<57r/3orais offset 
from this range by multiples of 2ir. The final phase a 



E{a) 313.4° 104.5° 0.0° 255.5° 46.6° 

Etc) 77.9° 20.6° 245.4° 20.6° 77.9° 

EJa) 240.0° 210.0° 300.0° 210.0° 240.0° 

Etb) 8.7° 94.3° 300.0° 325.7° 111.3° 

Etc) 111.3° 145.7° 300.0° 274.3° 8.7° 

TABLE I: Explicit pulse phases for composite pulses compris- 
ing five 7r pulses and depicted in Figs. El and El 



can then be varied in an attempt to minimize the second 
order error term, and it turns out that this term can in 
fact be entirely removed by choosing 



a = T2arcsin(^5/32). 



(43) 



Note that the two ± signs must be chosen with the op- 
posite signs as indicated above but the two solutions are 
fundamentally equivalent. The fidelity of this composite 
pulse in the absence of off-resonance errors is identical to 
that of Fi and BB1. However, as shown in Fig. Etc), the 
fidelity in the presence of off-resonance errors, although 
symmetric, is similar to that of Fi and much worse than 
that of BB1. Explicit pulse phases are listed in Table HI 



VII. OFF-RESONANCE ERRORS 

It is tempting to try to design an antisymmetric 
composite pulse which suppresses the second order off- 
resonance term in the same way as can be achieved for 
pulse strength errors, but this approach is not success- 
ful. The pulse phases needed to suppress the first order 
off-resonance error are <j)\ = arccos(ll/16) « 46.6° and 
02 = arccos(l/4) « 75.5°, with 03 = 0, 04 = —02 and 
05 = — 01 as usual. This pulse was briefly discussed 
by Odedra et al. who noted its poor performance [24] . 
Evaluating the propagator for this pulse shows that it 
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contains a second order error term, and (equivalently) 
the fidelity expression 

T=l-7r 2 f 4 /2 + 0(f) (44) 

is only correct to fourth order. 

The reason for this behaviour is not hard to find: the 
second order error term for off-resonance errors is more 
complex than for pulse strength errors. Expanding the 
single pulse propagator to second order gives 

V(tt, <j>) w 7r A 0+7r/ 2/i0 (45) 

with ji = 7r/ 2 /2. Applying the BCH relation in this case 
is fairly simple as the cross terms between these pulses 
are 0(/ 3 ), and so can be neglected, leaving the forms 

3 

as before, but 

A 2 = 4 £ E sin ^" - °* - ^ E ^ ( 47 ) 

where the left hand part of the expression is analogous to 
that for pulse strength errors, but the additional terms 
on the right hand side arise from the second order errors 
in the individual pulses. Designing a composite pulse 
with second-order tolerance of off-resonance errors re- 
quires suppression of all these terms, which will not be 
achieved by an antisymmetric composite pulse: the ad- 
ditional term depends on not (/>", and so is not sup- 
pressed along with the first order term. 

It is interesting to note, however, that the additional 
second order term takes the same form as the first order 
term for a pulse strength error. Thus any composite pulse 
with second order tolerance of off-resonance errors must 
have first order tolerance of pulse strength errors, that is 
the composite pulse must have simultaneous tolerance of 
both sorts of error. 



VIII. SIMULTANEOUS ERROR TOLERANCE 

To achieve simultaneous tolerance of pulse strength 
and off-resonance errors it is necessary to choose the com- 
posite pulse phases such that Eqns. [31] and 26] are both 
equal to zero. This might seem difficult to achieve, but is 
in fact straightforward as noted by Tycko et al. [9|, LLCj . 
The two sets of toggling frame phases, <pj and (j)" are al- 
ternately different by +7r/2 (for odd j) and — tt/2 (for 
even j). Thus if the phases are chosen such that Eq. [31] 
separately sums to zero for the odd terms and the even 
terms, then Eq. [46] will also sum to zero [24]. A little 
thought shows that this separate cancellation is necessary 
as well as sufficient for simultaneous error tolerance. 

This can, in fact, be achieved with a sequence of five 
pulses. The phase vectors describing the three odd- 
numbered pulses must form an equilateral triangle, while 
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0.0 0.5 1.0 1.5 a/7T 2.0 



FIG. 4: (Color online) Coefficient of the fourth order infi- 
delity terms for pulse strength error (F e ) and off-resonance 
error (Ff) in five pulse composite pulses with simultaneous 
compensation of both errors to first order. 

the vectors for the two even pulses must be antiparallel. 
Equivalently the three odd phases, <^_, (j)' 3 and <j)' b1 must 
differ by ±27r/3, providing two constraints on the pulse 
phases, while the even phases, <fi' 2 and <j>'^ must differ 
by 7r, providing an additional constraint. Finally the re- 
quirement $ = provides a fourth constraint, and so one 
of the five phases, say 02, can be varied at will. Solving 
the simultaneous equations gives the family of solutions 

((f)) = (7r + 2a, a, -tt/3, -5tt/3 - a, -7tt/3 - 2a) (48) 

where a can be chosen at will. This fa mily is well known 
in the context of inversion pulse s [p , [lQj , and has also 
been discussed as general rotors [24j, but its properties 
have not been fully explored. All members of the fam- 
ily suppress first order pulse strength and off-resonance 
errors, but the effect on second order errors depends on 
the value of a. 

Studies to date have largely concentrated on values of 
a which result in particularly simple pulse sequences. For 
example one can set = <p2 by choosing a = —tt. Off- 
setting all phase angles by 7r (reflecting the equivalence of 
7To and tt^ pulses), and rewriting everything in the range 
between and 2tt gives 

(</>) = (0, 0, 2tt/3, tt/3, 2tt/3) (49) 

which is the Si inversion pulse [H, [To|- Alternatively one 
can set <pi — </>5 by choosing a = — 57r/6 to get the time 
symmetric sequence 

(0) = (4tt/3, 7tt/6, 5tt/3, 7tt/6, 4tt/3). (50) 

Although this composite pulse may not appear familiar, 
it is in fact closely related to the Knill pulse I40l which is 
frequently used for dynamic decoupling [l2[ Il3|. While 
the original Knill pulse is frequently described as a 7To 
pulse followed by a z-rotation, it is better thought of as 
7r rotation around some other axis in the x?/-plane. 

A more interesting approach is to find the propagator 
fidelity of the family of pulses under both pulse strength 
and off-resonance errors. Expanding these to fourth order 
gives 

^ = l-^(f)V e + 0(e 6 ) (51) 




1.0 -1.0 -0.5 0.0 J 0.5 1.0 



FIG. 5: (Color online) Fidelity achieved by some composite pulses with n — 5 designed to suppress both pulse strength and 
off-resonance errors: (a) the symmetric Knill type pulse, (b) pulse with suppression of second order pulse strength errors, (c) 
pulse with suppression of second order off-resonance errors. The inmost contour is drawn at an infidelity of 10 -4 . 



and 



lf 4 F f 



O(f) 



(52) 



where the coefficients of the fourth order terms are 



1 1 =F 12 cos a + 4 cos 2a + 4 V3(sin 2a =p sin a) (53) 



and the =p symbols are minus signs for F e and plus signs 
for Ff. In fact these are essentially the same function, 
with a simply offset by tt. (Note that the different pre- 
factors in Eqns [51] and [52] reflects the different way in 
which the two errors are parameterized.) As the graphs 
in Fig. [4] show it is possible to have very significant con- 
trol over the size of these error coefficients. Interestingly 
the two pulses discussed above appear to use quite poor 
choices of a; in particular the Knill pulse occurs at a 
global maximum in the pulse strength error and a local 
maximum in the off-resonance error. However the Knill 
pulse does have the advantage of being time symmetric, 
and so its response to off-resonance errors is frequency 
symmetric, as shown in Fig. [5] 

There are obvious alternative approaches. First, con- 
sider the points where the two curves cross, which occur 
at a = 27r/3 and a = 5tt/3. At these points the sen- 
sitivities of the composite pulse to pulse strength and 
off-resonance errors are in some sense balanced. More 
interestingly, note that the two curves touch the x-axis, 
and so it is possible to find values of a such that one 
of the two coefficients is set to zero. This corresponds 
to finding a pulse sequence which eliminates the second 
order term in the corresponding error. The solution for 
removing second order pulse strength errors occurs at 



a = ± arccos[(3 =F V 13) /8] « 94.3° or 325.7° (54) 

while for removing second order off-resonance error the 
solution occurs at 



a = ±arccos[(-3TV13)/8] « 145.7° or 274.3°. (55) 



The two choices in each case give equivalent composite 
pulses; note that as expected from the discussion below 
Eq. [53] sequences optimised for pulse strength errors and 
off-resonance errors have values of a differing by 180°. 
Fidelity plots for these new composite pulses are shown 
in Fig. [5j confirming the expected behaviour, while some 
explicit pulse phases are listed in Table [I] 



These composite pulses also illustrate some general 
principles. It is possible to remove first and second or- 
der pulse strength errors without using an antisymmetric 
pulse: antisymmetry is a sufficient but not a necessary 
condition. This is important, as it is easy to show that 
no antisymmetric composite pulse with n = 5 can sup- 
press both pulse strength and off-resonance errors: anti- 
symmetric pulses must have <fi' 2 = (/> 4 , which contradicts 
the requirement that these toggling frame phases must 
be separated by tt. Once antisymmetry is abandoned 
it is possible to develop a composite pulse with n = 5 
which suppresses first and second order off-resonance er- 
rors, and this composite pulse also suppresses first order 
pulse strength errors. 



Finally I return to the different parametrization of 
pulse strength and off-resonance errors briefly discussed 
above. This is a general phenomenon, and as a result 
the inner contours in fidelity plots will be a factor of 
7r/2 wider in composite pulses optimized for off-resonance 
errors than in corresponding pulses optimized for pulse 
strength errors. Similarly the behaviour of the outer con- 
tours is dominated by the simple observation that, in the 
absence of off-resonance errors, the fidelity must always 
fall to zero at e = ±1, while the equivalent limits for 
off-resonance errors occur at / = ±y/3. For this reason 
fidelity plots for high order composite pulses will always 
look better in the off-resonance dimension than in the 
pulse strength dimension. 
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IX. SEQUENCES WITH SEVEN PULSES 

This approach can be extended to longer pulses, and 
remains fairly straightforward for the case n = 7. In 
this case there are three even pulses, which must form an 
equilateral triangle in the toggling frame, and four odd 
pulses which must form either a rhombus or a degener- 
ate arrowhead. The case of a rhombus leads to the two 
constraints = (/)[+ tt and <j)' 7 = 3 + 7r, while the case 
of an arrowhead requires 3 = (/)[+ tt and <j)' 7 = <j)' b + tt. 
Together with the two constraints on the even pulses and 
the usual constraint that <I> = this leaves two free pa- 
rameters which can be adjusted to fine-tune the pulse 
sequence. 

Solving for the case of the rhombus with two free pa- 
rameters leads to the phases 



(0) = (a, 2/?, 0, -27T/3, -tt/3 - 
-2tt/3 - 2a + 4(3, 



a - 



2/?, 
-2a 



3/3) (56) 



where a and f3 can be chosen at will. As before the 
fidelity can be written as a Taylor series in e and /, with 
fourth order coefficients that are functions of both a and 
f3. A pulse with simultaneous suppression of second order 
pulse strength and off-resonance errors can be found by 
solving for F e = Ff = 0, but as both coefficients are 
non-negative it is sufficient to solve the simpler equation 
F e + Ff = 0. Even this function is difficult to solve 
analytically, but it can be easily investigated numerically. 

Plotting F e + Ff as a function of a and f3 reveals 
four distinct minima where the function is equal to zero. 
These all occur at points where f3 = a — 27r/3, and im- 
posing this constraint makes the equation easy to solve. 
The solutions occur at 



. 1 U±Vu 

a = ± arccos ±-\ 

2 V 2 



(57) 



where the first two ± signs must take the same value, but 
the third sign can be varied independently. 

These four roots come in two pairs, which lead to 
equivalent pulse sequences, so there are only two gen- 
uinely distinct solutions. Of these two solutions one is 
much better than the other, as it has smaller higher or- 
der errors. The best result occurs at a ~ 192.8°, and its 
behaviour is shown in Fig.[6fa); the explicit pulse phases 
are given in Table HU These solutions were obtained as- 
suming that the even phases in the toggling frame in- 
crease around the equilateral triangle, but the results are 
entirely equivalent if decreasing angles are used instead. 

Next consider the case where the four odd pulses form 
an arrowhead in the toggling frame. Once again the equa- 
tions can be solved with two free parameters, but in this 
case it is not possible to achieve simultaneous suppres- 
sion of the second order pulse strength and off-resonance 
errors, and I do not pursue this possibility further. 

The behaviour of the symmetric Knill pulse in the case 
n = 5 suggests that there is some value in imposing over- 
all symmetry on the composite pulse sequence. This can 



01 02 


03 


04 


05 


06 


07 


Eta) 192.8° 145.7° 


72.8° 


240.0° 


252.8° 


145.7° 


12.8° 


Ep>) 252.5° 265.0° 


97.5° 


170.0° 


97.5° 


265.0° 


252.5° 


EJc) 72.5° 265.0° 


277.5° 


170.0° 


277.5° 


265.0° 


72.5° 



TABLE II: Explicit pulse phases for composite pulses com- 
prising seven ty pulses and depicted in Fig. [6] 



be achieved by imposing a single additional constraint, 
as symmetrising any pair of pulses is sufficient, together 
with the other constraints, to symmetrize the entire se- 
quence. Symmetric pulse sequences have antisymmetric 
phases in the toggling frame, and so 0' 5 = — 3 . This col- 
lapses the distinction between rhombus and arrowhead 
sequences, and the solutions can be found by requiring 
that 03 = 0i ± 7r. For the case of the plus sign the 
solution is 



(0) = (a, 2tt/3 + 2a, 7tt/3 + 3a, IOtt/3 + 4a, 

7?r/3 + 3a, 2tt/3 + 2a, a) (58) 

and the free parameter a can be varied in order to op- 
timize the suppression of either pulse strength or off- 
resonance errors. In particular, choosing 



(3 - >/6l) /16 



(59) 



completely removes the second order pulse strength error, 
while choosing 



a 



arccos 



(V61 -3) /16 



(60) 



completely removes the second order off-resonance error. 
The wider behaviour of these pulses is shown in Fig. [6f b) 
and (c), and explicit pulse phases are given in Table ITT1 
As usual these plots show better visible performance for 
off-resonance errors than for pulse strength errors, and 
values of a for pulse sequences optimised for off-resonance 
errors differ from those optimised for pulse strength er- 
rors by 180°. 

For the case of a minus sign, so that 3 = <\>\ — 7r, the 
general solution is 

(0) = (a, 2tt/3 + 2a, tt/3 + a, -2tt/3, 

tt/3 + a, 27r/3 + 2a, a) (61) 

As before it is possible to find sequences optimized for 
either pulse strength or off-resonance error, but in this 
case optimising for one sort of error seems to lead to 
large second order error terms of the other kind (they 
do, of course, continue to suppress all first order errors). 
The resulting composite pulses thus quite sharply favour 
one type of error over the other, and are not considered 
further here. 

Finally I consider what can be achieved with antisym- 
metric seven pulse sequences. As for n = 5 it is not pos- 
sible to simultaneously remove first order pulse strength 
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FIG. 6: (Color online) Fidelity achieved by various composite pulses with n — 7: (a) second order suppression of both pulse 
strength and off-resonance errors, (b) symmetric pulse optimized for pulse strength errors, (c) symmetric pulse optimized for 
off-resonance errors. The inmost contour is at an infidelity of 10 -5 . 



and off-resonance errors, as the symmetry of the toggling 
frame phases prevents this. In particular the requirement 
that 06 = 02 prevents the formation of an equilateral tri- 
angle from the even phases. 



X. SEQUENCES WITH NINE PULSES 

Tackling the general case with n = 9 is significantly 
more difficult than the lower numbers, as it is now dif- 
ficult to use geometric insight to make initial progress. 
The four even pulses must form a rhombus or arrowhead 
in the toggling frame, but the five odd pulses must form 
a pentagon, for which there are no simple geometric re- 
strictions. It is, however, possible to make some progress 
by imposing particular symmetries on the problem. 

The antisymmetric case was studied by Odedra et at. 
[24j . An antisymmetric pulse with n — 9 has only four 
controllable phases, and as the toggling frame phases 
must be symmetric there is no distinction between the 
rhombus and arrowhead cases. A solution, which they 
call ASBO-9, occurs when 

(0) (4a + 3a + 2^, 2a + a + tt, 0, 

— a — 7r, —2a — ijj, —3a — 2ip, —4a — (62) 

where ip = arccos(— 1/4) as usual, and a can be varied 
at will [24|. All such composite pulses suppress both 
first and second order error terms arising from both pulse 
strength and off-resonance errors, and higher order error 
terms can be partially controlled by the choice of a. 

It is not possible to completely remove either third or- 
der error term by varying a, but the third order pulse 
strength error can be minimized by choosing a ~ 308.0°, 
while the off-resonance term can be minimized by choos- 
ing a ~ 128.0°; these values of a differ by 180° as usual. 
Alternatively the two error terms can be balanced by 
choosing a = The first two solutions correspond to 





01 02 03 


04 


05 


07>6 


Ha) 


162.4° 162.4° 313.4° 


284.5° 


0.0° 


A 


13b) 


46.6° 255.5° 255.5° 


75.5° 


0.0° 


A 


13c) 


282.1° 339.5° 339.4° 


159.4° 


114.6° 


S 



TABLE III: Explicit pulse phases for some composite pulses 
comprising nine tt pulses and depicted in Fig. Q phases 0j>6 
can be obtained by noting that each pulse is either antisym- 
metric (A) or symmetric (S). 



the pulse sequences ASBO-9(Bi) and ASBO-9(^), while 
the latter two are the pulse sequences ASBO-9 (7A) and 
ASBO-9(7B) respectively [24]. The performance of the 
last two pulse sequences is shown in Fig. [7(a) and (b), 
and the phase angles are listed in Table IIII1 

Next I consider the design of symmetric pulse se- 
quences. Solving the constraint equations is quite com- 
plex in this case, but the solution is found to be 

(0) = (a, /3, /3, p - tt, 2/3 -2a, P- tt, f3, f3, a) (63) 

with 



P = 2a + arccos[-(l + 2cosa)/2] 



and 



s[(4 - Vl0)/4] 



(64) 



(65) 



The resulting phases are listed in Table IIII1 and the per- 
formance of this composite pulse is shown in Fig. [7(c). 



XI. OTHER APPROACHES 

It is useful to compare these results with other ap- 
proaches that have been used to develop composite pulses 
with simultaneous error tolerance. An early example is a 
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FIG. 7: (Color online) Fidelity achieved by various composite pulses with n = 9 giving simultaneous suppression of first and 
second order pulse strength and off-resonance errors: (a) ASBO-9(7A), (b) ASBO-9(7B), (c) a new symmetric pulse. The 
inmost contour is at an infidelity of 10 -6 . 



sequence of nine tt pulses described by Alway and Jones 
[22[ | with phases 

(0) = (0, 3^, 3^, ^, 7T - ^, tt + ^, (66) 

where -0 = arccos(— 1/4) as usual. This composite pulse 
suppresses first and second order error terms for pulse 
strength errors and first order off-resonance errors, but 
this can be achieved with a five pulse sequence as dis- 
cussed above. It should, however, be noted that this 
pulse was developed as an illustration of a general ap- 
proach to designing composite pulses, in which individ- 
ual errors are removed one by one [20-22], rather than 
being optimized for these particular properties. 

A more interesting approach is the family of concate- 
nated composite pules developed by Bando, Ichikawa, 
Kondo and Nakahara [30|, IfUj , which seek to achieve si- 
multaneous error correction by nesting correction of off- 
resonance errors within correction of pulse strength er- 
rors, or vice versa. Once again simultaneous error toler- 
ance can be achieved, but the results are no better than 
for sequences of five tt pulses. However these concate- 
nated pulse can be extended to pulses with rotation an- 
gles other than n, and here they may be more useful. A 
similar approach has been explored by Merrill and Brown 

0. 

If, however, simultaneous error tolerance is not re- 
quired then previous approaches to composite pulse de- 
sign are likely to be preferable in many cases. In particu- 
lar for the pulse strength errors the fully time symmetric 
BB1 sequence suppresses first and second order errors 
with a total sequence length equivalent to five tt pulses, 
and at no cost in increased sensitivity to off-resonance 
errors. If even more effective error suppression is re- 
quired then symmetrised W n pulses [19] can be used. 
Similarly, for off-resonance errors the CORPSE pulse se- 
quence [5|,|42[ can be used to suppress the first order error 
and almost (but not entirely) remove the second order er- 
ror, with a total sequence length equivalent to between 



four and five tt pulses, and at no cost in increased sen- 
sitivity to pulse strength errors. All of these composite 
pulses have the advantage that they can be generalized 
to other rotation angles. 



XII. CONCLUSIONS 

It is clear that it is possible to design quite short com- 
posite not gates with excellent simultaneous tolerance 
of pulse strength and off-resonance errors: comparing 
Fig. 0(c) with Fig. [2](a) shows that the sequence of nine 
tt pulses has a similar region of parameter space inside 
the 10 -6 contour as the naive single tt pulse has inside 
the contour at 10 -2 . 

The insight from NMR studies that it is useful to con- 
centrate on the geometric form of the error, rather than 
proceeding blindly with algebraic minimisation, is cer- 
tainly correct for pulses up to n = 9. Beyond this it 
is less clear how such geometric constraints should be 
applied, as suitable closed equilateral polygons can be 
formed in many ways. However it is less obvious that 
the second NMR insight, that antisymmetric composite 
pulses which suppress first order errors also automati- 
cally suppress second order errors, is as useful, and for 
many purposes in QIP symmetric composite pulses, de- 
signed by combining geometric and algebraic elimination 
of error terms, are likely to be more appropriate. 

One further advantage of antisymmetric pulses de- 
signed to tackle pulse strength errors is that they can 
be iteratively nested [19] to produce composite pulses 
with exceptionally broad error tolerance. This approach 
also works with pulses designed for simultaneous error 
tolerance; in particular ASBO-9 composite pulses can be 
nested to produce a composite pulse which removes both 
error terms up to ninth order, so that the infidelity is 
eighteenth order in both e and /. This behaviour can ul- 
timately be traced to the fact that these antisymmetric 
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pulses can suppress multiple even order error terms [7j. Acknowledgments 
However, nested ASBO-9 pulses comprise 81 separate tt 
pulses, and such long composite pulses are not considered 

here. I thank Ron Daniel for helpful conversations. 
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